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Abstract 

The  present  study  is  focused  on  the  optimization  of  the  proton  exchange  membrane  fuel  cell  (PEMFC).  A  simplified  conjugate-gradient  method 
(SCGM)  is  employed  to  seek  the  optimal  combination  of  the  design  parameters,  including  the  channel  width  ratio  (A),  the  porosity  of  GDL  (sgdl)> 
and  the  porosity  of  the  catalyst  layer  (eCat)-  In  this  report,  the  concept  of  the  method  is  described  in  detail  and  the  results  of  the  optimization  process 
are  discussed.  The  results  show  that  the  present  approach  is  capable  of  providing  an  optimal  design  under  specific  operation  conditions. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

In  the  past  several  decades,  fossil  fuels  have  caused  environ¬ 
mental  pollution  problems,  such  as  the  greenhouse  effect  and 
acid  rain.  Due  to  the  increasing  rate  of  consumption  of  global 
petroleum,  finding  other  more  efficient  and  cleaner  substitutive 
power  generation  systems  has  become  an  issue  of  great  concern. 
Among  the  possible  alternative  power  generation  systems,  pro¬ 
ton  exchange  membrane  (PEM)  fuel  cells  have  received  much 
attention  from  energy-related  researchers  and  industries  world¬ 
wide.  Since  in  PEM  fuel  cells  electric  power  is  generated  by 
directly  transforming  the  chemical  fuel  energy  into  electricity 
in  electrochemical  reactions,  these  devices  exhibit  a  higher  effi¬ 
ciency  compared  to  traditional  power  generators.  Also,  without 
ignition  and  combustion  during  electricity  generation,  PEM  fuel 
cells  are  free  of  severe  pollution  and  noise.  In  addition,  PEM 
fuel  cells  do  not  have  the  flaws  of  corrosion  of  the  mechani¬ 
cal  elements  and  material  problems  associated  with  the  working 
medium.  Recently,  PEM  fuel  cell  technology  has  been  developed 
more  actively  than  ever  in  the  automobile  and  3C  (computer, 
communication,  consumer)  electronics  industries.  The  com¬ 
mercialisation  of  PEM  fuel  cells  into  mass  production  can  be 
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accelerated.  However,  more  efforts  are  still  needed  to  explore 
the  physical  characteristic  features  as  well  as  to  improve  the 
performance  of  PEM  fuel  cells  in  the  coming  years. 

A  number  of  reports  have  been  presented  in  the  past  several 
decades;  however,  only  those  closely  related  to  the  present  study 
are  mentioned  in  the  following. 

Hental  et  al.  [  1  ]  considered  the  effects  of  carbon  plate  calcina¬ 
tions  at  high  temperature  on  the  PEMFC  efficiency.  The  authors 
claimed  that,  with  a  flow  channel  width  of  1  mm  and  the  thick¬ 
ness  of  the  collector  ribs  of  around  2  mm,  the  performance  of 
their  fuel  cell  was  obviously  elevated. 

Hontanon  et  al.  [2]  used  CFD  software  to  study  the  effects 
of  GDL  permeability  and  the  ratio  of  the  flow  channel  to  the  rib 
widths  on  the  PEMFC  performance.  The  results  showed  that  a 
decrease  in  the  permeability  of  the  gas  diffusion  layer  (GDL) 
reduces  cell  performance.  It  was  also  found  that  the  influence  of 
the  GDL  permeability  on  the  cell  performance  is  stronger  than 
that  of  the  flow  channel  width  ratio. 

Yoon  et  al.  [3]  similarly  studied  the  influences  of  the  flow 
channel  width  on  the  PEMFC  efficiency.  Their  research  showed 
that  the  gas  diffusion  effect  plays  a  more  important  role  than  the 
electric  conductivity  effect.  As  the  flow  channel  width  is  kept 
constant,  the  performance  of  fuel  cell  is  higher  with  a  smaller 
rib  width. 

Gurau  et  al.  [4]  developed  a  two-dimensional  numerical 
model  to  study  the  changes  in  the  gas  concentration  and 
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Nomenclature 

D 

mass  diffusivity  (m2  s-1) 

F 

Faraday  constant  (96,487  C  mol-1) 

I 

current  density  (Am-2) 

j 

transfer  current  density  (Am-3) 

J 

objective  function 

ib 

base  width  of  the  flow  distributor  (m) 

k 

gas  channel  width  (m) 

k 

shoulder  width  (m) 

Mo2 

molar  mass  of  oxygen  (32.0kgkmol-1) 

n 

iteration  step 

nQ 

number  of  electrons 

P 

power  output  per  unit  area  (W  m-2) 

R 

universal  gas  constant  (8.314  J  (mol  K)-1) 

AjDL 

thickness  of  gas  diffusion  layer  (m) 

tC  at 

thickness  of  catalyst  layer  (m) 

T 

temperature  (K) 

V 

potential 

Greek  letters 

a 

transfer  coefficient  for  the  reaction 

p 

step  size 

r 

conjugate-gradient  coefficient 

£Cat 

porosity  of  catalyst  layer 

£GDL 

porosity  of  gas  diffusion  layer 

T] 

over-potential  (V) 

K 

reaction  rate  constant 

X 

channel  width  ratio,  /c//b 

p 

density  of  oxygen  (kgm-3) 

a 

electric  conductivity  (Q~l  m-1) 

^Cat 

tortuosity  of  catalyst  layer 

Jgdl 

tortuosity  of  gas  diffusion  layer 

0 

electric  potential  (V) 

CO 

oxygen  mass  fraction 

(L>0 

oxygen  mass  fraction  at  inlet 

electrochemical  reaction  rate  by  including  the  electrochemical 
kinetics  of  the  catalyst  layer  and  the  convection-diffusion  effect 
in  the  GDL.  Results  show  that  as  the  porosity  of  the  GDL  is 
increased,  the  entry  velocity  is  elevated.  They  also  found  that 
the  cell  performance  increases  with  the  cell  temperature. 

Wang  and  Savinell  [5]  studied  the  effects  of  porosity,  thick¬ 
ness,  Pt-loading,  CO  poisoning  of  catalyst  layer  on  the  cell 
performance.  In  this  report,  it  was  observed  that  the  catalyst 
layer  porosity  of  an  efficient  fuel  cell  is  between  0.4  and  0.7. 

Yan  et  al.  [6]  studied  the  influence  of  flow  channel  geom¬ 
etry  and  GDL  porosity  on  gas  transport.  Results  provided  in 
this  report  showed  that  a  larger  flow  channel  width  ratio  or  a 
higher  GDL  porosity  leads  to  a  better  cell  performance.  They 
also  pointed  out  that  an  increase  in  the  number  of  flow  channels 
results  in  a  more  uniform  current  density  distribution  and  thus  a 
better  cell  performance. 

Zhou  and  Liu  [7]  proposed  a  three-dimensional  fuel  cell 
model  to  study  the  influence  of  CO  toxins  and  the  variation 


of  hydrogen  concentration  in  anode.  Authors  claimed  that  some 
phenomena,  which  were  not  found  by  using  the  one-  and  two- 
dimensional  models,  are  worthy  of  investigation  by  the  three- 
dimensional  model.  They  also  carried  out  a  parametric  analysis 
of  the  effects  of  the  entry  velocity,  the  porosity  of  the  GDL,  and 
the  channel  width  ratio  on  the  cell  performance.  Based  on  this 
report,  it  was  observed  that  the  channel  width  ratio  and  the  GDL 
porosity  both  have  their  own  optimal  values. 

Kumar  and  Reddy  [8]  studied  the  improvement  of  the  perfor¬ 
mance  of  the  PEMFC  based  on  optimization  of  the  dimensions 
and  shape  of  the  bipolar  plates.  It  was  found  that  the  optimal 
dimensions  for  the  channel  width,  the  land  width  and  the  chan¬ 
nel  depth  were  close  to  1.5,  0.5  and  1.5  mm,  respectively. 

Hsu’s  work  [9]  was  concerned  with  the  analysis  of  the  trans¬ 
port  phenomenon  and  the  parametric  optimization  of  a  PEM 
fuel  cell.  This  study  was  focused  on  the  influence  of  conductiv¬ 
ity,  channel  width  ratio,  and  porosity  and  thickness  of  the  GDL. 
A  two-dimensional  half-cell  model  is  utilized  in  simulation  of 
the  transport  phenomenon. 

Ying  et  al.  [10]  investigated  the  effects  of  channel  configu¬ 
ration  on  an  air-breathing  fuel  cell’s  performance.  The  authors 
found  that  the  best  performance  can  be  obtained  in  the  cell  with 
a  3  mm  wide  cathode  channel  and  an  open  ratio  of  75.9%. 

Lum  and  McGuirk  [11]  carried  out  a  series  of  parametric 
studies  of  the  effects  of  variation  of  the  electrode  thickness, 
shoulder  width,  degree  of  permeability  and  oxidant  concentra¬ 
tion,  so  as  to  provide  a  clearer  understanding  of  how  changes  in 
these  parameters  affect  the  cell  performance. 

Berning  and  Djilali  [12]  presented  the  results  of  a  parametric 
study  based  on  a  three-dimensional,  non-isothermal  model  of  a 
PEM  fuel  cell.  The  authors  evaluated  the  influence  of  the  channel 
width  ratio  and  the  thickness  and  porosity  of  the  GDL. 

From  the  above  survey  of  related  literature,  it  is  noted  that 
most  of  the  studies  were  focused  on  numerical  simulation  of 
the  transport  phenomenon  in  a  PEM  fuel  cell,  or  at  most  on  the 
parametric  study  of  the  effects  of  physical  variables.  Only  a  few 
reports  dealt  with  the  optimization  of  a  set  of  design  parameters. 

In  practice,  the  task  of  optimization  searches  for  a  combi¬ 
nation  of  the  most  appropriate  values  of  a  number  of  design 
parameters  with  a  minimized  objective  function.  Nowadays, 
several  useful  optimization  methods  are  available,  for  exam¬ 
ple,  the  genetic  algorithm  [13],  the  conjugate-gradient  method 
(CGM)  [14],  the  Newton-Raphson  method  [15],  the  adaptive 
weighting  input  estimation  method  [16],  etc.  These  methods  are 
all  applicable  to  minimize  the  prescribed  objective  functions. 
Among  them,  the  CGM  method,  which  stems  from  perturbation 
principles,  transforms  the  optimal  problem  into  three  separate 
problems,  namely,  the  direct  problem,  the  sensitivity  problem, 
and  the  adjoining  problem.  The  solutions  of  the  three  problems 
give  the  updated  conjugate-gradient  direction  as  well  as  the  step 
size  towards  convergence.  Cheng  and  Chang  [17,18]  used  the 
concept  of  the  conjugate-gradient  method  to  develop  an  inverse 
method  for  predicting  the  internal  temperature  distribution  of 
PEM  fuel  cells.  The  results  show  that  the  outer  surface  tem¬ 
perature  data  can  be  used  to  predict  the  internal  temperature 
distribution  at  the  interface  between  the  carbon  plate  and  the 
MEA. 
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Under  these  circumstances,  in  the  present  study,  a  multi- 
parametric  optimization  for  the  PEMFC  design  is  investigated. 
Among  the  possible  influential  parameters,  the  channel  width 
ratio,  and  the  porosities  of  the  GDL  and  the  catalyst  layer  are  rec¬ 
ognized  as  the  influential  factors  that  may  affect  the  behaviour  of 
proton  exchange  and  electron  conduction;  therefore,  these  three 
parameters  are  selected  as  the  design  parameters  to  be  optimized 
in  this  work.  The  simplified  conjugate-gradient  method  (SCGM) 
proposed  by  Cheng  and  Chang  [14]  is  employed  herein  as  the 
optimization  search  scheme.  This  study  is  aimed  at  optimiza¬ 
tion  of  these  three  parameters  so  that  the  best  performance  of 
the  PEMFC  under  some  specified  geometry  and  operating  con¬ 
ditions  can  be  obtained.  However,  it  is  noted  that  the  present 
optimization  approach  is  not  limited  to  the  present  group  of 
designed  parameters.  Where  necessary,  more  designed  parame¬ 
ters  may  be  readily  applied. 

2.  Optimization  method 

The  present  optimization  method  consists  of  two  major 
numerical  solutions:  (1)  the  direct  problem  solver  and  (2)  the 
simplified  conjugate-gradient  method  (SCGM).  The  two  solvers 
are  briefly  described  as  follows. 

2.7.  Direct  problem  solver 


1 .  The  fuel  cell  operates  steadily. 

2.  The  total  pressure  and  temperature  are  uniform  and  constant 
in  the  solution  domain. 

3.  The  gases  are  assumed  to  be  compressible  ideal  gases,  fully 
saturated  with  water  vapor. 

4.  The  oxygen  flux  is  zero  at  the  interface  between  the  catalyst 
layer  and  the  membrane. 

5.  The  gas  transport  from  the  flow  channel  to  the  catalyst  layer 
is  only  by  diffusion. 

6.  The  gas  diffusion  layer  and  the  catalyst  layer  are  considered 
to  be  isotropic  porous  media. 

7.  Owing  to  high  ohmic  resistance  of  the  membrane,  the  inter¬ 
face  between  the  catalyst  layer  and  the  membrane  is  assumed 
to  be  insulated  against  electrons. 

The  channel  width  ratio  (A)  is  defined  by  a  ratio  of  the  channel 

width  /c  to  the  carbon  plate  width  4>  as 


2.7.7.  Catalyst  layer 

In  the  catalyst  layer,  the  oxygen  is  consumed  due  to  the 
electro-chemical  reaction  therefore,  a  two-dimensional  trans¬ 
port  equation  for  the  oxygen  mass  fraction  (cocat)  in  the  catalyst 
layer  can  be  expressed  as 


Fig.  1  shows  a  schematic  of  a  two-dimensional  half-cell 
model  for  the  cathode  side  employed  for  simulation  of  the 
transport  phenomena  involved.  From  top  to  bottom,  the  model 
consists  of  the  catalyst  layer,  the  gas  diffusion  layer,  the  flow 
channel,  and  the  carbon  plate  which  is  developed  with  the  fol¬ 
lowing  simplifications  and  assumptions: 


BC6 


b 


i 

i 

tCat 

Catalyst  Layer 

r 

BC5 

J 

i 

Igdi. 

Gas  Diffusion  Layer 

BC7 

BC3 

BC4 

i  i 

i 

tch 

Flow  Channel 

BC8 

b 

BC2 

y 

y 

r 

- ► - 

0  x  BC1 


BC9 


Fig.  1.  Schematic  of  the  two-dimensional  simulation  domain  of  half-cell 
PEMFC  (cathode). 


^Cat 

Cat 


where  po2  stands  for  the  oxygen  density,  Do2  the  oxygen  diffu- 
sivity,  F  the  Faraday  constant,  j  the  transfer  current  density,  nQ 
the  number  of  electrons  in  transfer  current  density  relation,  and 
£Cat  and  Teat  represent  the  porosity  and  the  tortuosity  factor  of 
the  catalyst  layer,  respectively. 

The  cathode  transfer  current  density  can  be  expressed  as 


j  —  UqFkcoqxp 


where  k  is  the  reaction  rate  constant,  a  the  transfer  coefficient  for 
the  reaction,  rj  the  over-potential,  R  the  universal  gas  constant, 
and  T  is  the  cell  temperature. 

Due  to  the  electro-chemical  reactions  activated  in  the  catalyst 
layer,  there  exists  a  source  term  in  the  Poisson  equation  for  the 
electric  potential.  That  is 


a  (  eft  &t> Cat\ 

Tx\ac*^r) 


(eff  #Cat  \ 

l  Cat  dy  ) 


=  (1  -  £Cat)TCa,OCat 


(5) 

(6) 


where  creat  denotes  the  electric  conductivity  of  the  catalyst  layer 
and  0cat  the  electric  potential  of  the  catalyst  layer. 
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2.7.2.  Gas  diffusion  layer  ( GDL ) 

In  the  GDL,  only  the  oxygen  diffusion  takes  place  without 
chemical  reactions.  The  two-dimensional  transport  equation  for 
the  oxygen  mass  fraction  &>gdl  in  the  GDL  can  be  expressed 
as 

9  ( ,,  neff  9"GDL  \  9  (  neff  9"GDL  \  „ 

dx  V  °2 °GDL ”17“ )  +  dy  V  °2 °GDL )-°  (7) 

with 


D 


eff 

GDL 


_  orGDL 

—  £gdl 


where  £gdl  and  tgdl  are  the  porosity  and  the  tortuosity  factor 
of  the  GDL,  respectively. 

Since  in  the  GDL  there  is  no  electrochemical  reaction  taking 
place,  the  electric  potential  equation  possesses  no  electric  source 
term  and  is  expressed  as 


d 

dx 


^GDL 

aGDL 


#GDL  \ 
dx  ) 


^GDL 

aGDL 


90GDL \ 

) 


OqDL  =  (1  “  ^GDL ) TCat  G'GDL 


(9) 

(10) 


where  <tgdl  is  the  electric  conductivity  and  0gdl  is  the  electric 
potential  in  the  GDL. 


2.1.3.  Carbon  plate 

It  is  essential  to  investigate  the  effects  of  gas  channel  width 
on  the  transport  phenomena  of  the  cathode  reactants  and  on  the 
ohmic  resistance  of  the  whole  module. 

The  carbon  plate  serves  as  a  conductor  for  the  electric  cur¬ 
rent.  The  electric  potential  equation  in  the  carbon  plate  can  be 
expressed  as 


9  (  #P\  ,  9 

dx  Lp  dx  )  +  dy 


(11) 


where  <rp  is  the  electric  conductivity  and  </>p  is  the  electric  poten¬ 
tial  in  the  carbon  plate. 


2.1.4.  Boundary  conditions 

The  physical  conditions  at  the  boundaries  numbered  B 1  to 
B9  in  Fig.  1  are  described  as  follows: 

The  boundary  conditions  for  oxygen  mass  fraction  ( co )  are 


daJG  DL 

dy 

9&>gdl 

dx 


at  boundaries  BC4; 
at  boundaries  BC7  and  BC9 


The  boundary  conditions  for  electric  potential  (0)  are: 


Catalyst  layer: 
90Cat 


dy 


=  0  at  boundary  BC6; 


0Cat  =  0GDL  at  boundary  BC5 ; 

0Cat  _  q  at]30un(tariesBC7  andBC9 

dx 

GDL: 

^eff  l  ^GDL  =  ^Cat  ^  boundary  BC5 . 


dy 


dy 


90GDL 

dy 


=  0  at  boundary  BC3 : 


0GDL  =  0p  at  boundary  BC4; 

90GDL 

- =0  at  boundaries  BC7  and  BC9 

dx 

Carbon  plate: 

0p  =  0  at  boundary  BC 1 ; 


df 


P  9y 
#P 


p  =  <Tqdl  ^>GDL  at  boundary  BC4; 


dy 


dx 

9</>P 

dy 


=  0  at  boundaries  BC7,  BC8,  andBC9: 


=  0  at  boundary  BC2 


2.1.5.  Operation  potential 

With  the  solutions  for  the  electric  potential  distribution 
yielded  using  the  above  model,  the  cell  potential  (V)  can  be 
calculated  by  the  following  equation: 

V  =  Voc  ~  0  -  A0  (12) 


Catalyst  layer: 


9&>Cat 

dy 


=  0  at  boundary  BC6; 


coCai  =  &>gdl  at  boundary  BC5; 


9&>Cat 

dx 


=  0  at  boundaries  BC7  and  BC9 


GDL: 


D 


eff  9<£>GDL  ^eff  9<^Cat 


GDL 


=  D 


dy  Cat  dy 

(Ugdl  =  ojq  at  boundary  B  3 


at  boundary  BC5; 


where  Voc  denotes  the  ideal  cell  potential,  r\  is  for  the  over¬ 
potential  of  the  cathode,  and 


p  (06  -  <pi)dx 
L 


(13) 


denoting  the  average  ohmic  loss  in  terms  of  the  difference 
between  the  electric  potentials  at  boundaries  B6  and  Bl, 
respectively.  During  the  optimization,  the  cell  voltage  V  and 
cell  current  I  are  varied  accompanying  iterative  changes  in  the 
design  parameters. 

It  is  noted  that  the  present  two-dimensional  model  is  rather 
straightforward  and  not  able  to  include  the  effects  of  heat  and  liq¬ 
uid  water  accumulation.  A  more  powerful  model  will  certainly 
be  required.  However,  the  simplicity  of  the  model  is  acceptable, 
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given  that  only  the  channel  width  ratio  and  the  porosities  of  the 
GDL  and  the  catalyst  layer  are  taken  into  account  and  the  opti¬ 
mization  is  emphasized  rather  than  the  complexity  of  the  model. 


2.2.  Simplified  conjugate- gradient  method  (SCGM) 

In  this  study,  the  simplified  conjugate-gradient  method 
(SCGM)  [10]  is  employed  to  optimize  the  values  of  the  design 
parameters  in  search  of  optimal  performance.  The  objective 
function  in  conjunction  with  the  optimization  process  is  defined 
in  the  following: 

J(Ak)=j,  k  =  1,2,3  (14) 

where  (k=  1,2,  and  3)  represent  the  channel  width  ratio  (A), 
the  porosity  of  the  GDL  (£gdl),  and  the  porosity  of  the  catalyst 
layer(£cat)«  The  cell  power  output  (P)  is  expected  to  reach  a  max¬ 
imum  as  the  magnitude  of  the  objective  function  /  is  minimized. 
In  the  SCGM  method,  the  values  of  the  step  size  fik  during  the 
searching  procedure  are  fixed.  In  this  manner,  the  complexities 
in  the  mathematical  manipulation  and  the  constraint  of  the  form 
of  the  objective  function  can  be  alleviated.  The  initial  guess  for 
each  design  variable  is  made  first,  and  in  the  successive  steps  the 
conjugate-gradient  coefficients  and  the  searching  directions  are 
evaluated.  The  new  design  variables  are  continuously  improved. 
When  the  objective  function  reaches  a  minimum,  the  goal  of  the 
optimization  is  reached  and  the  searching  procedure  is  then  ter¬ 
minated. 

The  procedure  for  the  iterative  optimization  process  with  the 
SCGM  method  is  presented  as  follows: 


(1)  Make  the  initial  guess  for  the  designed  variables  A&,  assign 
the  values  to  the  step  sizes  /3&. 

(2)  Specify  all  boundary  conditions,  and  then  solve  the  direct 
problem. 

(3)  Calculate  the  objective  function  J(Ak,  k  =  1-3).  If  the  conver¬ 
gence  criterion  is  satisfied  as  the  objective  function  reaches 
a  minimum  in  the  iterating  process,  the  solution  process  is 
terminated.  Otherwise,  proceed  to  step  (4). 

(4)  Calculate  the  gradient  functions  of  the  function  /: 


dJn  _  a  r 

dAk  A Ak  ’ 


k  =  1-3 


(15) 


(5)  Calculate  the  conjugate-gradient  coefficients  yj!  (k=  1-3): 


(dJ/dAkf 

(dJ/dAk)n~l 


k=  1-3 


(16) 


For  the  first  step  with  n  =  0,  y)  =  0. 

(6)  Calculate  the  searching  directions  ^  (k=  1-3): 


tn+ 1 


k  =  1-3 


(17) 


(7)  Update  new  design  variables  by 

Ar'=Al-fk%»,  k=  1-3  (18) 


Table  1 

Variables  used  in  the  present  analysis 


Parameters 

Notation 

Value 

Temperature  (K) 

T 

333 

Reaction  rate  constant  (mol  (m3  s)-1) 

K 

1.5  x  10“5 

Oxygen  mass  fraction  at  inlet 

(i>0 

0.17 

Gas  channel  width  (m) 

lc 

1.0  X  10“3 

Channel  shoulder  width  (m) 

h 

1.0  x  10“3 

Thickness  of  gas  diffusion  layer  (m) 

tG  DL 

5  x  10“4 

Thickness  of  catalyst  layer  (m) 

*Cat 

5  x  10“5 

Porosity  of  GDL 

£gdl 

0.5 

Tortuosity  of  GDL 

^GDL 

1.5 

Porosity  of  catalyst  layer 

£Cat 

0.3 

Tortuosity  of  catalyst  layer 

^Cat 

1.5 

Laraday  constant  (C  mol-1) 

F 

96487 

Universal  gas  constant  (J  (molK)-1) 

R 

8.314 

Molar  mass  of  oxygen  (kg  kmol-1) 

mo2 

32.0 

Number  of  electrons 

ne 

2 

Transfer  coefficient  for  reaction 

a 

1 

Electric  conductivity  of  catalyst  layer  (£2-1  m-1) 

tfCat 

190 

Electric  conductivity  of  GDL  (£2-1  m-1) 

O'GDL 

537 

Electric  conductivity  in  carbon  plate  (£2-1  m-1) 

4000 

£2-1m-1,  <7p  =  4000  £2-1  m-1,  ^gdl  =  5  x  10  4  m,  teat  =  5  x  10-5  m,  Z  =  0.5, 
£gdl  =  0-5,  and  scat  =  0.3.  (a)  Polarization  curve  and  (b)  power  curve  of  a  PEM 
fuel  cell. 


and  then  return  to  step  (2). 
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3.  Results  and  discussion 

3.1.  Preliminary  study 

In  the  computation,  the  activated  over-potential  acts  as  an 
independent  variable  to  replace  the  external  loads.  By  doing  so, 
the  polarization  diagram  of  the  fuel  cell  can  be  determined.  The 
power  density  of  the  fuel  cell  is  the  product  of  the  potential  V 
and  the  current  density  I  as 

P=VI  (19) 

The  model  described  in  the  preceding  section  is  solved  by  a 
finite  volume  method  with  a  uniformly  distributed  grid  having 
101  and  511  nodes  in  the  x-  and  y-  directions,  respectively.  The 
values  of  the  geometrical  and  physical  variables  considered  in 
the  present  study  are  listed  in  Table  1 . 

Fig.  2  shows  the  polarization  and  the  power  curves 
for  a  typical  case  at  the  conditions  of  cr  cat  =  190  £2_1  m_1, 


c)rGDL  =  537  £2  1  m  1 ,  crp  =  4000  £2  1  m  1,^gdl  =  5x10  4  m, 
teat  =  5  x  10-5  m,  a  =  0.5,  £qdl  -  0-5,  and  £cat  =  0.3.  For  a  lower 
over-potential,  the  weaker  electrochemical  reaction  leads  to 
a  lower  power  density.  However,  when  the  over-potential  is 
increased,  the  reaction  is  stronger  so  that  the  power  density 
reaches  a  maximum  at  I-  7400  A  m2. 

Sensitivity  analysis  of  the  three  optimized  design  param¬ 
eters  with  respect  to  the  cell  voltage  has  been  carried 
out  and  Fig.  3(a)-(c)  are  plotted  to  show  the  results.  The 
cases  shown  in  Fig.  3(a)-(c)  are  at  treat  =  190  m-1, 

ctgdl  =  537  £2_1  m-1,  crp  =  4000  £2_1  m-1,  *gdl  =  5  x  10-4  m, 
Aat-5  x  10_5m.  Fig.  3(a)  shows  the  effects  of  the  channel 
width  ratio  (A)  for  £gdl  =  0.5  and  £cat  =  0.3,  Fig.  3(b)  shows 
the  effects  of  porosity  of  GDL  (£gdl)  for  A  =  0.5  and  £cat  =  0.3, 
and  Fig.  3(c)  conveys  the  effects  of  porosity  of  catalyst  layer 
(scat)  f°r  A  =  0.5  and  £gdl  =  0.5.  It  is  found  that  in  the  ranges 
of  parameters  considered,  the  cell  voltage  is  more  sensitive  to  a 
change  in  the  channel  width  ratio  or  the  porosity  of  the  GDL. 


Fig.  3.  Dependence  of  cell  voltage  on  designed  parameters  at  crCat  =  190  m_1,  (Tgdl  =  537  m-1,  crp  =  4000  m-1,  foDL  =  5  x  10-4  m,  feat  =  5  x  10-5  m. 

(a)  Effects  of  channel  width  ratio  (A.)  for  £gdl  =  0.5  and  £cat  =  0.3.  (b)  Effects  of  porosity  of  GDL  (£gdl)  for  A  =  0.5  and  £cat  =  0.3.  (c)  Effects  of  porosity  of  catalyst 
layer  (£cat)  f°r  A  =  0.5  and  £gdl  =  0.5. 
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However,  it  is  relatively  less  sensitive  to  the  porosity  of  the  cat¬ 
alyst  layer. 

Fig.  4  further  shows  the  effects  of  the  channel  width  ratio 
(A)  on  the  performance  of  a  PEM  fuel  cell  at  r\-  0.25  V, 
?GDL  =  5  X  10-4  m,  ?Cat  =  5  X  10-5  m,  eGDL  =  0.5,  and  eGat  =  0.3. 
It  is  observed  that  indeed  the  channel  width  ratio  has  a  subtle 
influence  on  the  cell  performance.  Larger  channel  width  ratios 
may  enhance  the  fuel  gas  transport  as  well  as  the  chemical  reac¬ 
tion  in  the  catalyst  layer;  however,  a  larger  channel  width  leads 
to  smaller  ribs  and  hence  results  in  an  increase  in  the  ohmic 
impedance.  With  these  two  opposing  effects,  there  must  be  an 
optimal  value  for  the  channel  width  ratio  (A),  at  which  the  perfor¬ 
mance  reaches  a  peak  value.  In  Fig.  4,  the  results  show  that  the 
largest  power  of  6572  W  m-2  is  obtained  at  the  optimal  channel 
width  ratio  of  A  =  0.65  under  the  conditions  considered  in  this 
case. 

Fig.  5  shows  the  effects  of  the  GDL  porosity  (£gdl)  on  the 
performance  of  a  PEM  fuel  cell  at  rj  =  0.25  V.  A  higher  GDL 
porosity  implies  a  higher  ohmic  impedance  and  results  in  a 
reduction  in  the  electric  current.  However,  on  the  other  hand,  a 
higher  GDL  porosity  also  implies  a  stronger  diffusion  of  reactant 
gas  in  the  GDL.  Fig.  5  shows  that  the  optimal  value  of  £gdl  is  0.6 
at  rj-  0.25  V,  £cat  =  0.3,  ^gdl  =  5  x  10_4m,  tcat  =  5  x  10_5m, 
and  A  =  0.5.  At  the  optimal  value  of  £gdl,  the  power  density 
approaches  6843  W  m~2. 

The  effects  of  the  catalyst  layer  porosity  (£cat)  on  the 
performance  of  a  PEM  fuel  cell  at  r\-  0.25  V,  £gdl  =  0.5, 
^gdl  =  5  x  10-4  m,  teat  =  5  x  10-5  m,  and  A  =  0.5  are  plotted  in 
Fig.  6.  Except  for  the  case  at  very  low  or  very  high  porosity 
£Cat,  the  effects  of  the  catalyst  layer  porosity  are  not  as  obvi¬ 
ous  as  those  of  the  GDL  porosity.  This  is  attributed  to  the  fact 
that  the  catalyst  layer  is  relatively  thin,  and  hence  the  effects 
of  £cat  on  the  performance  are  limited.  The  variation  of  the 
power  output  is  only  less  than  1  %  as  the  porosity  £cat  is  in  the 
range  of  0.3-0. 8.  Even  so,  the  optimal  porosity  of  the  catalyst 
layer  £cat  is  also  obtained.  For  the  case  at  rj  =  0.25  V,  £gdl  =  0.5, 


Fig.  4.  Effects  of  channel  width  ratio  k  on  performance  of  a  PEM  fuel  cell  at 
7]  =  0.25  V,  tGDL  =  5  x  10~4  m,  tc at  =  5  x  10-5  m,  £Gdl  =  0.5,  and  £Cat  =  0.3. 


Fig.  5.  Effects  of  GDL  porosity  £gdl  on  performance  of  a  PEM  fuel  cell  at 
r]  =  0.25  V,  £cat  =  0.3,  foDL  =  5  x  10-4  m,  tc  at  =  5  x  10-5  m,  and  k  =  0.5. 

tGDL  =  5  x  10-4  m,  feat  =  5  x  10-5  m,  and  A  =  0.5,  the  optimal 
£Cat  is  0.5  and  the  corresponding  power  of  6099  W  m-2  can  be 
observed. 

3.2.  Optimization  of  PEMFC 

In  the  parametric  studies  considered  above,  only  one  param¬ 
eter  is  changed  while  others  are  fixed.  In  fact,  multi-parameter 
optimization  is  more  meaningful  from  the  design  viewpoint. 
Fig.  7  shows  the  trajectory  of  the  multi-parameter  optimization 
process  approaching  the  optimal  design.  The  optimal  combi¬ 
nation  of  A,  £gdl,  and  £cat  are  sought  in  the  parameter  space 
(A,  £gdl,  £Cat)«  The  optimization  process  is  started  with  the 
initial  guess  of  A  =  0.65,  £gdl  =  0.5,  and  £cat  =  0.3.  It  can  be 
observed  in  Fig.  7  that  the  SCGM  method  effectively  brings  the 
searching  process  to  approach  an  optimal  state.  In  this  case, 


Fig.  6.  Effects  of  catalyst  layer  porosity  £Gat  on  performance  of  a  PEM  fuel  cell 
at  i]  =  0.25  V,  £Gdl  =  0.5,  fGDL  =  5  x  10~4  m,  fGat  =  5  x  10~5  m,  and  k  =  0.5. 
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Fig.  7.  Trajectory  of  the  optimization  process  for  a  PEM  fuel  cell  seeking  an 
optimal  set  of  designed  parameters  X,  sqdl,  and  £Gat- 


Initial  guess  1  Initial  guess  2 


Fig.  9.  Trajectories  of  the  optimization  processes  from  different  sets  of  initial 
guesses  for  X,  £Gdl,  and  £Gat- 


it  is  found  that  the  iteration  converges  after  about  400  iter¬ 
ations.  The  computation  time  required  is  approximately  80  h 
on  a  personal  computer  with  an  AMD-3GMHz  CPU.  Based 
on  the  results  shown  in  Fig.  7,  it  is  found  that  at  A.  =  0.55, 
£gdl  =  0.5913  and  £cat  =  0-307,  the  optimization  performance 
can  be  generated  with  a  maximum  power  of  6866.53  W  m~1 2 3 4  at 
0  =  0.25  V,  foDL  =  5  x  10-4  m,  feat  =  5  x  10~5 6  m.  The  variation 
of  the  power  density  during  the  optimization  process  is  shown 
in  Fig.  8.  Note  that  the  power  density  associated  with  the  initial 
design  is  approximately  2500  Wm~2.  The  power  density  then 
approaches  its  maximum  value  when  the  optimal  parameters  are 
obtained. 

To  check  the  robustness  of  the  searching  method  used  in  the 
present  work,  Fig.  9  shows  the  convergence  of  the  optimization 
process  from  three  different  initial  guesses  for  the  three  param¬ 
eters  A,  £gdl,  and  scat*  at  r)  =  0.25  V.  It  can  be  observed  in  Fig.  9 


that  the  present  optimization  method  leads  to  the  same  optimal 
set  of  parameters  with  different  initial  guesses. 

4.  Conclusions 

In  this  study,  an  optimization  analysis  employing  the  SCGM 
method  [10]  combined  with  a  two-dimensional  half-cell  model 
for  the  transport  phenomena  has  been  successfully  developed. 
The  optimization  of  design  parameters  A,  £gdl,  £cat>  has  been 
performed.  The  results  show  that,  at  the  conditions  of  A  =  0.55, 
£gdl  =  0.5913  and  £cat  =  0.307,  the  optimal  cell  performance  in 
terms  of  power  is  P  =  6866.53  W  m-2.  Extension  of  the  present 
method  to  one  with  more  parameters  is  a  worthwhile  work  for 
the  future,  and  a  more  comprehensive  model  would  also  be  inter¬ 
esting.  It  is  noted  that  the  present  optimization  approach  is  not 
limited  to  the  present  group  of  designed  parameters.  Where  nec¬ 
essary,  more  designed  parameters  may  be  readily  applied. 


Fig.  8.  Variation  in  power  output  during  the  optimization  process. 
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